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We systematically analyze the coherence length in even-even nuclei. The pairing coherence length 
in the spin-singlet channel for the effective density dependent delta (DDD) and Gaussian interaction 
is estimated. We consider in our calculations bound states as well as narrow resonances. It turns out 
that the pairing gaps given by the DDD interaction are similar to those of the Gaussian potential 
if one renormalizes the radial width to the nuclear radius. The correlations induced by the pairing 
interaction have in all considered cases a long range character inside the nucleus and decrease towards 
the surface. The mean coherence length is larger than the geometrical radius for light nuclei and 
approaches this value for heavy nuclei. The effect of the temperature and states in continuum is 
investigated. Strong shell effects are evidenced, especially for protons. We generalize this concept 
to quartets by considering similar relations, but between proton and neutron pairs. The quartet 
coherence length has a similar shape, but with larger values on the nuclear surface. We evidence the 
important role of proton-neutron correlations by estimating the so-called alpha coherence length, 
which takes into account the overlap with the proton-neutron part of the a-particle wave function. 
It turns out that it does not depend on the nuclear size and has a value comparable to the free a- 
particle radius. We have shown that pairing correlations are mainly concentrated inside the nucleus, 
while quarteting correlations are connected to the nuclear surface. 
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I. INTRODUCTION 

The concept of coherence has a general character being 
connected to the linear superposition of quantum states. 
Two-body coherence properties in nuclear structure are 
directly connected to the properties of low-lying collec¬ 
tive states. Collective excitations are microscopically de¬ 
scribed by a superposition of creation pair operators act¬ 
ing on the ground state, described by a coherent state 
within the Random Phase Approximation (RPA). The 
coherent state in this context is defined as an exponential 
excitation of products between pair operators acting on 
the vacuum state Q ■ It is well known that ground state 
properties of even-even nuclei are well reproduced by 
the pairing interaction M- The wave function within 
the Bardeen-Cooper-Schrieffer (BCS) pairing approach is 
also of a coherent type, i.e. an exponential excitation of 
the pair creation operators acting on the vacuum state. 

The spatial distribution of the two-particle density 
is very important in understanding nuclear correlations 
[§, ,f|. In particular, in Ref. [7] it was analyzed the rela¬ 
tion between coherence and chaotic properties of the nu¬ 
clear pairing. The coherence property is characterized by 
the so-called coherence length, defined as the root mean 
square distance averaged over the density. For superfluid 
nuclei this average is usually performed over the pairing 
density. In Refs. nm it was shown that this quantity is 
relatively large, comparable to the nuclear size inside the 


nucleus and decreases beyond the nuclear surface. This 
picture of an extended di-nuclear cluster can be under¬ 
stood in terms of the Pauli blocking, hindering nucleons 
to cluster together inside the nucleus and, therefore, the 
cluster loses binding and becomes larger. It is in contrast 
to the a-clustering phenomenon, which_ takes place in a 
narrow region close to the surface area 0 m , being con¬ 
nected to the very large binding energy of an a-particle 
moving in a low density region [13| . Thus, we expect that 
the corresponding correlation length estimated between 
proton and neutron pairs will have significantly smaller 
value. 

The finiteness of nuclear systems also has important 
consequences as far as thermal properties are concerned. 
Pairing correlations in finite nuclei do not vanish at some 
critical temperature, but they slowly decrease over sev¬ 
eral MeV 0, [jj| . This can be theoretically obtaine d by 
projecting the particle number in the BCS theory [16| . 
However, hints about such behaviour can be extracted in 
the unprojected BCS approach from the spatial proper¬ 
ties of the correlations. 

In this paper we will perform a systematic analysis of 
the pairing coherence length and a comparison to the sim¬ 
ilar quantity defined for quartets. In Section II we give 
the necessary theoretical background concerning pair¬ 
ing equations containing resonant states and coherence 
length. In Section III we perform a systematic analysis 
of the coherence length and in the last Section we draw 


2 


Conclusions. 


II. THEORETICAL BACKGROUND 
II. 1. Pairing equations 

In order to investigate two-body correlations we ex¬ 
pand the wave function of N + 2 particles in terms of the 
wave function of N particles as follows 

|^JV+2) = k\V N ) • (2.1) 


the operator k, connecting N + 2 with N systems, is 
called pairing density operator. In this case the expan¬ 
sion coefficient 


X e = ~(BCS n+2 I [at ® al] 0 \BCS N ) 
%/2 j + 1 


(2.3) 


is given in terms of BCS amplitudes as follows 


= u ( jv + 2 ) w ( jv ) [. 
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(2.4) 


We will consider in our calculations the spherical approx¬ 
imation. Thus, the operator a\ creates a single particle 
(sp) eigenstate of the spherical mean field potential with 
standard quantum numbers e = (eZj). In the configura¬ 
tion representation one has 


(r, s|a| JO) = r, s) = + e (r) ® xi (s) 

+<J r ) = ip e (r)i l Yi M (f) = ^p-i l Yi^{f) , (2.2) 


where yjr) is the radial wave function and the rest of 
the notation is standard. 

The operator k in Eq. m is called within the decay 
theory two-particle formation amplitude. In the absence 
of two-body correlations, when the wave functions are 
Slater determinants, this relation is nothing else than 
the Laplace expansion of the (N+2) x (N+2) normalized 
determinant in terms of NxN times 2x2 normalized de¬ 
terminants. 

The most important two-body correlation beyond the 
mean field in even-even nuclei is given by the pairing in¬ 
teraction. We will describe such systems within the stan¬ 
dard BCS approach, where the averaged particle number 
is conserved, separately for protons and neutrons. Thus, 
both wave functions in Eq. m have a BCS ansatz and 
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We will consider in our basis bound sp states with neg¬ 
ative energy, as well as relatively narrow sp resonances 
with positive energy. Relatively narrow resonances are 
similar to bound states and can be normalized to unity 
in the internal region, but at large distances they behave 
like outgoing waves 


<Pe(r) 


M e — 


(+) 


M s Mi G,(r)+.'F,(r) ^ 


in terms of spherical Hankel functions for neutrons and 
Coulomb-Hankel functions for protons. The coefficients 
M e are called scattering amplitudes and their squared 
values are proportional to sp partial decay widths. 

The states in continuum play an important role on 
pairing correlations, especially for nuclei close to the drip 
lines js, Il7i42ll|. For nuclear structure calculations the 
background contribution is not relevant and only rela¬ 
tively narrow resonant states are important [22j, [231 • A 
very good approximation for BCS calculations is to ne¬ 
glect the finite resonance width, i.e. to treat the reso¬ 
nances as bound-like states [24[. We label bound states 
by a and resonances with positive energy by r. We treat 
proton and neutron pairing separately; for a given isospin 
index the generalized system of BCS equations for gap 
parameters A Q , A r and number of particles N is 
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( 2 . 6 ) 


where A is the chemical potential and the potential ma¬ 
trix elements V a ,p are computed according to Eq. (2.6) 
of [25|. 

We will investigate pairing in excited nuclei by us¬ 
ing the temperature-dependent equations with anoma¬ 


lous and normal densities, respectively 


\a e ae) = u e v e tanh— 

(■ a t a e ) = vl + {u\ - u e 2 )/(e' 3Be + 1) . 


(2.7) 
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II.2. Pairing coherence length 


The coherence length is defined as follows 


The two-body operator entering the pairing density 
m can be written in the configuration representation. 
By using the recoupling from j-j to the L-S scheme, one 
obtains spin-singlet and spin-triplet components. Our 
calculations have shown that the largest contribution is 
given by the spin-singlet component, given the following 
expression 


«(ri, r 2 ) = E [^e(ri) ® ^(r 2 )] 0 

Z ' r i T'n 


( 2 . 8 ) 


in terms of two-particle azimuthal harmonics 
yi (cos 9) = [i l Yi{h)®i l Yi(r 2 )\ 0 

V21 + T 


47T 


-Pi (cos 9) , 


(2.9) 


where 9 is the angle between particle radii and the ex¬ 
pansion coefficient is given by 


= x e \/j + i((ZZ) 0, 


1 1 
2 2 


0;0| Z- R Z- j;0) 


( 2 . 10 ) 

in terms of LS-jj recoupling brackets. By expanding the 
sp wave function with respect to the harmonic oscillator 
(ho) basis 

<M r ) 


</>nil( r ) = 4>nl M*%(f) , 


dp) 


( 2 . 11 ) 


where [3 = M^u/H is the standard ho parameter, and 
by using the Talmi-Moshinski transformation to rela¬ 
tive r = ri — r 2 and center of mass (c.o.m.) coordinate 
R = (i’i + r 2 )/2 one obtains the following expansion 


n(r,R,9) = ^2 f\(r,R)y x (cos9) , 

A 

with expansion coefficients given by 

fx(r,R) = E^4a /2) W ®nx(R) 


( 2 . 12 ) 

(2.13) 


nN 


where 


QnN\ — ^ ^ ( c nie c n2li (nXN\, 0|niZn 2 Z, 0) . 

e n 1 n 2 

(2.14) 

Here, the bracket denotes the standard Talmi-Moshinsky 
recoupling coefficient. By averaging over the angle 9 we 
get 

1 f 1 

R 2 (r,R) = — / K 2 (r,R,9)dcos9 

2 J -1 
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1 7(2)(1?) 
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= fj dr r 2 w(r , R) , (2.16) 

in terms of the integrals 

POO 

/ ( Ri?) = / dr r 2p H 2 (r, R) 

= E 


X NN' 

x y yn.N\y n 'N'\ 


^^ QnNxQn 


dr r 2p (t^l 2) {r)(j) ( P/ 2 \r) . 


(2.17) 


Let us finally mention that the quantity x e , defined by 
Eq. (12.41) . is also called ’’anomalous” density, while the 
quantity 


Ve = v 2 


(2.18) 


is called ’’normal” density. Therefore k. defined by Eq. 
(12.81) . can be called ’’anomalous” coherence length, while 
a similar quantity kq defined by using the ’’normal” den¬ 
sity is called ’’normal” coherence length. 


II.3. Quarteting correlations 

We will investigate quarteting correlations in medium 
and heavy a-decaying nuclei, where the valence protons 
and neutrons occupy different major shells. The standard 
assumption to build a quartet from two protons and two 
neutrons in such nuclei is to consider proton and neutron 
pairing separately j26[ [23] . Therefore the system of N n + 
2,1V,, + 2 nucleons can be expressed in terms of N n . N v 
nucleons in a factorized way as follows 

14/iv jr -t- 2 ,iv 1/ -t- 2 } = k- K k v \'l’ n^n„) , (2-19) 

where n T is defined by Eq. GED- Thus, the quartet 
wave function is a product between proton and neutron 
two-body wave functions (12.121) . Anyway, calculations 
in infinite nuclear matter suggest that a-clusters can oc¬ 
cur only at relative small nuclear densities compared to 
the equilibrium value and the proton-neutron correla¬ 
tions play an important role [13|. Thus, an a-particle 
can be formed only in the surface region where the nu¬ 
clear density diminishes and proton-neutron correlations 
become relevant. This situation can be simulated by a 
proper modification of the single particle mean field by 
adding a gaussian interaction in the surface region [28} 
and still by keeping the factorized ansatz (12.191) . This 
can expain why an a-particle can be formed from two 
protons and two neutrons lying in different major shells. 
This additional ansatz of the single particle mean field 
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was recently confirmed by microscopic calculations [29j 
and fission-like theory [301 ]. Anyway, this modification 
is important in order to reproduce the absolute value of 
the half-life, but has a minor influence on the coherence 
length. 

In order to describe quartets we introduce the relative 
and c.o.m. coordinates for proton, neutron and proton- 
neutron systems, respectively: 

„ r±+r 2 

r?r = ri - r 2 , Rtt =-^— 

I? r 3 + r 4 

r„ = r 3 - r 4 , R„ = —-— 

Fa = Rtt R^ , Ra = ^^±1^ , (2.20) 

where we labeled by 1, 2 proton and by 3, 4 neutron 
coodinates. The internal a-particle wave function is 
given by the product between the lowest proton, neutron 
and proton-neutron ho orbitals 

V’a = ^ /2 \rM^ 2 \rM^\r a ) , (2.21) 

where /3 a ss 0.5 fm~ 2 is the free a-particle ho parameter 
measured by electron scattering experiments [261 ]. This 
parameter is about 2-3 times larger than the similar sp 
ho parameter in heavy a-emitters, due to the fact that 
a-particle is a very bound object. 

We will describe quarteting correlations between pro¬ 
ton and neutron pairs by overlapping the relative coordi¬ 
nates to the corresponding components of the a-particle 
wave function (12.211) . We will proceed in two steps. 

A. Quarteting correlation length 

Let us first consider only the overlap with respect to 
proton and neutron relative coordinates tv, tv, by keep¬ 
ing free the internal proton-neutron coordinate r a . Thus, 
we consider independent from each other proton and 
neutron pairs by neglecting proton-neutron correlations. 
Therefore we can define the quarteting density in anal¬ 
ogy to the pairing density, but between the proton and 
neutron pairs (instead of fermions): 

n q {R v , R^) = (K w (ri,r 2 )|^oo“ /2) ( r 7r)} 

x (^(r 3 ,r 4 )|4o“ /2) (^)) • (2-22) 

By recoupling the product between proton and neutron 
pairs (12.121) to the relative and c.o.m. pair coordinates 
one obtains for the leading monopole component the fol¬ 
lowing relation 

K q (r a ,R a )tt K W(r a ,R a )= ^ G n (N n ) G V (N V ) 

N n ,N u 

x ^(ncO-ZVaO; 0\N % ON V 0; 0)</>jfJ 0 ( r a )<j>^Q(Ra) , 

n a 

(2.23) 


in terms of Moshinsky brackets and the proton/neutron 
monopole G-coefficients (12.141) 

Gr(N T ) = ^S^NM&\rr)\<t>^ a/a \r T )) 

n T 

T = 7T, V . (2.24) 

It does not depend on angles and therefore one can de¬ 
fine the quarteting coherence length £ q (R a ) without any 
additional angular average (12.151) by using in Eqs. (12.161) 
and (12.171) the quarteting density squared K q (r a , R a ). 

B. Alpha coherence length 

The next step is to consider proton-neutron correla¬ 
tions. They are described by the corresponding part in 
the a-particle wave function (12.211) given by </>qq a \r a )- 
In order to account for the narrow proton-neutron spa¬ 
tial distribution in the free a-particle one defines the so- 
called alpha coherence length £ a (R a ) by using the alpha 
density 

Ra )^00 (^"a) ; (2.25) 

in performing the integrals (12.171) . 

Let us finally mention that the integral of the alpha 
density over the relative proton-neutron coordinate 

n OO 

T{R a )= / n a {r a ,R a )r 2 a dr a , (2.26) 

Jo 

defines the formation amplitude and its square describes 
the probability to find an a-particle in the quartet wave 
function [l^, [26|. 



FIG. 1. Paring gap defined by the first two lines of Eq. (ESI) 
versus e in 48 Cr for DDD potential (squares) and Gaussian 
potentials with ro= 2 fm (circles) and ro = Rn (diamonds). 


III. NUMERICAL APPLICATION 

We analyzed all even-even nuclei with 20 < Z < 100 
and known experimental pairing gaps, determined by the 
binding energies of neighbouring nuclei [(jit . 
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TABLE I. Proton quantum numbers, sp spectrum, decay widths and gap parameters for the Gaussian, renormalized Gaussian 
and DDD interactions in 48 Cr, given by the diagonalization of the Woods-Saxon mean field with universal parametrisation 

m- 


No. 

l 

2 j 

e (MeV) 

T (MeV) 

A 2fm (MeV) 

A 4 . 5 fm(MeV) 

Addd(MbV) 

1 

0 

1 

-28.911 

- 

3.114 

1.354 

0.724 

2 

1 

3 

-20.837 

- 

3.173 

1.810 

1.482 

3 

1 

1 

-18.638 

- 

3.121 

1.739 

1.436 

4 

2 

5 

-12.118 

- 

2.908 

2.131 

2.387 

5 

0 

1 

-8.349 

- 

2.454 

1.795 

1.728 

6 

2 

3 

-7.488 

- 

2.886 

2.047 

2.351 

7 

3 

7 

-3.079 

- 

.261 

2.224 

2.246 

8 

1 

3 

0.322 

0.000 

1.349 

1.356 

1.076 

9 

1 

1 

2.403 

0.046 

1.149 

1.133 

0.962 

10 

3 

5 

4.101 

0.024 

2.114 

2.003 

2.139 

11 

4 

9 

5.874 

0.055 

1.389 

1.893 

0.996 


For the nuclear mean field we used a standard Woods- 
Saxon potential with universal parametrization fEE . We 
considered in our sp basis all bound states and resonances 
in continuum up to e mQX = 10 MeV with a sp decay width 
r < 1 MeV. As an example we give in Table [Q the proton 
sp spectrum for 48 Cr. Here, there are given level number, 
angular momentum, twice the total spin, sp energy, decay 
width of sp states in continuum and pairing gaps for the 
interactions considered below. 



FIG. 2. Proton coherence length divided by geometrical 
radius versus c.o.m. radius in 48 Cr computed with ’’normal” 
(a) and ’’anomalous” densities (b) for DDD potential (solid 
line) and Gaussian potentials with ro= 2 fm (long dashes) 
and r o = Rn (short dashes). 


We solved the BCS equations (12.611 separately for pro¬ 
tons and neutrons with two widely used types of nucleon- 
nucleon pairing interactions: 

I. Gaussian interaction 


It is defined by the following ansatz: 

v(ri 2 ) = -v 0 e~^ ri2 ^ r ^ 2 , (3.1) 

depending on the relative radius 772 . Here, the width pa¬ 
rameter ro=2 fm corresponds to the spin-singlet ’’bare” 
value in the free space. The corresponding value of the 
effective potential strength vo is determined by the gap 
parameter at the Fermi level, which should be equal to 
the experimental value. 

II. Density dependent delta (DDD) interaction 


It is known that the strength of the effective pairing 
interaction depends upon the local density fl7l El , given 
by the following phenomenological ansatz [19( 


u(r, r') = uqS(t — r') < 1 — X 


Pn( r) 


-| 7 > 


( 0 ) 

Pn 


(3.2) 


in terms of the nuclear density pN- The value X=1 cor¬ 
responds to the surface DDD interaction. 

As an example, in Fig. (T] we plotted the pairing gap 
(12.61) versus sp energy for 48 Cr, given in Table Q] Here, 
circles correspond to the Gaussian interaction in the free 
space with ro=2 fm. Notice large values for states below 
the Fermi level. The gaps given by DDD interaction with 
X = 7=1 are plotted by squares and the values below the 
Fermi level are significantly smaller that the Fermi gap. 

It is interesting to point out that a very similar be¬ 
haviour has the Gaussian interaction where the width 
parameter is renormalized to the geometrical nuclear ra¬ 
dius (in fm) ro = R n = 1.2A 1 / 3 . A di-nuclear cluster 
inside nuclear matter has different properties with re¬ 
spect to the free space. It considerably loses the binding 
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FIG. 3. The integrand of the proton coherence length versus 
the relative radius radius in 48 Cr, computed with ’’normal” 
(a) and ’’anomalous” densities (b) for different c.o.m. radii. 
Here, we used the Gaussian interaction withn ro=2 fm. 


property due to the Pauli blocking, becoming larger and 
therefore the effective pairing interaction has a more ex¬ 
tended shape. Above the Fermi sea we obtained similar 
values in all cases. 

In Fig. [2] (a) we plotted the proton coherence length 
given by Eq. (12.161) divided by the nuclear radius Rn, 
as a function of the ratio between c.o.m. and nuclear 
radius R/Rn in 48 Cr. Here we used the ’’normal” den¬ 
sity while in Fig. [2] (b) we used the ’’anomalous” density. 
Notice that that all cases, plotted by different symbols 
explained in caption, have very similar shapes. Thus, the 
coherence length is not sensitive to the radial shape of the 
interaction. The ’’normal” coherence length is equal to 
the nuclear radius in the internal region and diminishes 
by a factor 0.5 on the surface. The ”anomalous” coher¬ 
ence length has a similar shape, but with twice larger 
internal value. This picture is very different from the 
dependence of the two-body wave function versus c.o.m. 
radius, which is peaked on the nuclear surface (25]. 

In order to better understand the behaviour of the co¬ 
herence length we plotted in Fig. [3] (a) the integrand of 
the ’’normal” correlation length wq( r,R), given by the 
second line of Eq. (12.161) . versus the relative radius r for 
three values of the c.o.m. radius R= 4.5 fm (solid line), 
2 fm (long dashes) and 0 fm (short dashes). Here, we 
used the ” bare” version of the Gaussian interaction. No¬ 
tice that the three curves have a similar shape, strongly 
peaked around 2 fm. We obtain completely different plots 
for the integrand of the ’’anomalous” coherence length 
w(r,R). They are given in Fig. [3] (b). The distribution 
corresponding to the c.o.m. radius on surface f?=4.5 fm 
(solid line) is peaked around the free singlet value of the 
Gaussian width i.e. r =2 fm. On the contrary, the dis¬ 
tribution corresponding to a smaller radius R= 2.5 (long 
dashes) is peaked around a much larger value r=7 fm. 

Our conclusions are in agreement with Ref. ;§(, where 
in Fig. 5 the ’’anomalous” coherence length of the pairing 
interaction was estimated within the more sophisticated 
Hartree-Fock-Bogoljubov (HFB) approach, by using the 


FIG. 4. (a) Proton coherence length versus c.o.m. radius 

for different chemical potentials A=-2.96 MeV (solid line), 0 
MeV (long dashes) and 0.96 MeV (short dashes) in 48 Cr. 

(b) Neutron coherence length versus c.o.m. radius for dif¬ 
ferent chemical potentials A=-13.75 MeV (solid line), 0 MeV 
(long dashes) and 0.88 MeV (short dashes) in 48 Cr. 



N N 


FIG. 5. Strength parameter of the Gaussian interaction, 
corresponding to ro=2 fm, versus neutron number for proton 
(a) and neutron systems (b). 


Gogny force for Ni isotopes. The shape is similar, pre¬ 
dicting a mean coherence length of about 6 fm in the 
internal region and decreasing as one approaches the nu¬ 
clear surface and reaching the value of 2 fm just outside 
the nucleus. 

Most of the exotic nuclei close to the drip lines have the 
last nucleon in continuum. Therefore we investigated the 
dependence of the coherence length on the Fermi level, by 
changing the real part of the Woods-Saxon potential. We 
plotted in Fig. 0] (a) the proton coherence length versus 
c.o.m. radius in 48 Cr, for different values of the chemical 
potential. One sees that it increases by increasing the 
chemical potential. This effect is stronger for neutrons, 
as seen in Fig. [U(b), due to the absence of the Coulomb 
barrier. Therefore, in exotic nuclei close to drip lines the 
nucleons become more correlated. 

Then we performed a systematic analysis of the 
’’anomalous” coherence length (by simply calling it co¬ 
herence length) for even-even nuclei with 20 < Z < 100. 
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FIG. 6 . Ratio (£)/Rn, corresponding to a Gaussian interac- FIG. 8 . Logarithm of the ratio (£)/Rn versus logarithm of 

tion with ro =2 fm, versus neutron number for proton (a) and the mass number for protons (a) and neutrons (b). 

neutron systems (b). 



FIG. 7. Ratio (£) /Rn, corresponding to the DDD interaction 
(IT21) with X = 7 = 1, versus neutron number for proton (a) 
and neutron systems (b). 


In Fig. [5] (a) we plotted the effective strength vq as a 
function of the neutron number for protons correspond¬ 
ing to Gaussian interaction with ro=2 fm. The isotope 
chains are connected by solid lines and magic numbers 
are indicated by vertical lines. Different regions are 
plotted by open squares (20 < Z < 28), filled squares 
(28 < Z < 50), open circles (50 < Z < 82) and filled 
circles (82 < Z < 100). As a general trend we remark a 
strong decreasing behaviour with the increase of the neu¬ 
tron number. We notice a remarkable feature, namely it 
has almost the singlet ” bare” value in the free space Vo ~ 
35 MeV for very light nuclei. The strength strongly de¬ 
creases up to wo ~ 20 MeV for heavy nuclei, except the 
regions around magic numbers. In Fig. [5] (b) we give a 
similar plot for neutrons. Notice that in this case shell 
effects are stronger. 

In Fig. [G] (a) we analyzed the mean coherence length 
(£) for protons, corresponding to the Gaussian interac¬ 
tion with the free value of the width parameter ro =2 fm 
as a function of neutrons. The ratio of this quantity to 
the nuclear radius decreases from 1.4 for light nuclei up 
to around unity for heavy nuclei. In Fig. [G] (b) we give 
similar results for neutrons. As a general trend, the co¬ 
herence length is larger for neutrons due to the absence 
of the Coulomb barrier, but the shell effects are stronger 



R (fm) 


FIG. 9. Proton coherence length versus c.o.m. radius in 
220 Ra computed for T=0 (solid line), T=0.5675 MeV < T c 
(long dashes), for a Gaussian potential with ro= 2 fm. 


for protons. 

We then investigated the density dependent pairing in¬ 
teraction given by Eq. (13.21) with A" = q=l in Fig. [7] 
It turns out that the ratio {Q/Rn has similar gross fea¬ 
tures, but with more pronounced shell oscillations. The 
fact that the coherence length for neutrons is larger is 
confirmed. It is interesting to notice the linear correla¬ 
tion between logio(t;) and log^A , plotted in Fig. [ 8 ] for 
the Gaussian pairing interaction with r$= 2 fm. 

In order to investigate the behaviour of the coherence 
length for excited states, in Fig. [9] we analyzed the role of 
the temperature. Firstly, we give for 220 Ra the coherence 
length versus the pair c.o.m. radius for T = 0 and just 
below the ’critical’ temperature T c ss 0.57 MeV (here the 
gap decreases below 10 -3 MeV). The pairing coherence 
length shows very little change in shape up to T c . The 


































































































FIG. 10. (a) Quarteting coherence length in 220 Ra versus 

c.o.m. radius. 

(b) Same as in (a) but for alpha coherence length. 
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FIG. 11. Averaged alpha coherence length versus mass num¬ 
ber. 


strongest variation appears in the internal region, while 
on the surface, where the pairs are strongly coupled Q, 
there is indeed almost no change. As a measure of the 
pairing correlations, the coherence lenght would appear 
to indicate a gradual transition to the normal state with 
increasing temperature, as its behaviour is similar to that 
of the pairing gap in the particle number conserving case 

mi). 

Our purpose is to compare the pairing and quarteting 
coherence lengths. First we analyzed the quarteting co¬ 
herence length, by using the quarteting density (12.221) , for 
the a-emitter 220 Ra as a function of c.o.m. radius in Fig. 
E0( a). One notices a similar qualitative behaviour com¬ 
pared to the pairing coherence length, but the absolute 
values are larger on the nuclear surface. Our calcula¬ 
tions have shown that the temperature practically does 
not change this dependence. 






FIG. 12. The two terms I < ' P \R), p=2 (solid line) p=l 
(dashed line) given by Eq,. (12.171) . defining the pairing co¬ 
herence length for protons (a), neutrons (b), quarteting co¬ 
herence length (c) and alpha coherence length (d) versus the 
c.o.m. radius. 


It turns out that the proton-neutron correlations, given 
by the overlap with the corresponding proton-neutron 
part of the a-particle wave function (12.251) . completely 
change this picture. One sees from Fig. [TO] (b), where 
we plotted the alpha correlation length versus c.o.m. ra¬ 
dius, that the values oscillate around the value of the 
geometrical radius of the a-particle R a . Thus, our anal¬ 
ysis confirm the crucial role played by proton-neutron 
correlations in the formation of the a-particle. Finally, 
in Fig. [FT] we plotted the mean value of the alpha co¬ 
herence length for even-even a-emitters above A = 100. 
It has a quasi-constant value around 1.7 frn. Small local 
maxima correspond to regions above double magic nuclei 
132 Sn and 208 Pb. 

In order to better understand the difference between 
pairing and quarteting correlations we plotted in Fig. [12] 
the two terms I^ p \ p=2 (solid line) and p=l (dashed line) 
given by Eq. (12.171) versus the c.o.m. radius. The two 
terms reach their maximal values for the pairing case (left 
panels) at R=0, while for the quarteting case (right pan¬ 
els) the maxima are centered around the surface region. 
The pairing coherence length for protons (a) and neu¬ 
trons (b) is given by the ratio between solid and dashed 
curves which obviously decreases with increasing c.o.m. 
radius. Quarteting coherence length is given by the ra¬ 
tio between solid and dashed lines in Fig. Q7] (c) which 
have slighly shifted broad maxima located below the nu¬ 
clear surface. Although the two terms have completely 
different shapes compared to the pairing case, their ratio 
plotted in Fig. [TO] (a) is also a decreasing function with 
respect to the c.o.m. radius. 
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The alpha coherence length, given by the ratio of the 
two curves in Fig. [12] (d), deserves special attention. 
These curves have very narrow maxima centered at the 
same point on the nuclear surface. Moreover, it turns out 
that the two curves are almost proportional and therefore 
their ratio leads to the quasiconstant value in Fig. fTOlfb). 
close the a-particle geometrical radius = 1.2 4 1 * / 3 4 5 6 * 8 9 10 11 ss 
1.9 fin. Notice that the shape of the curves in Fig. [T2] 
(d), peaked on the nuclear surface, is similar to the stan¬ 
dard a-particle formation probability, given by the inte¬ 
gral (12.261) squared fl2l |28[. 


IV. CONCLUSIONS 

In conclusion, we have performed in this paper a sys¬ 
tematic analysis of the pairing coherence length in the 
spin-singlet channel for various types of pairing interac¬ 
tion. We compared the DDD potential to the Gaussian 
interaction. We considered in our calculations bound 
states as well as narrow resonances. 

As a very important conclusion we have shown that, 
by considering the singlet ’’bare” value of the width pa¬ 
rameter tq=2 fm, the strength parameter reproducing 
the gap parameter for light nuclei is close to the singlet 
value in the free space vq ~ 35 MeV and decreases up to 
vq ~ 20 MeV for heavy nuclei. We have shown that the 
’’renormalized” Gaussian interaction with a larger width 
parameter than its free value ro=2 fm (equal to the nu¬ 
clear radius) has similar properties to the commonly used 
density dependent pairing potential. 

It turns out that the pairing coherence length has 
similar properties for all considered interactions. It is 


[1] P. Ring and P. Schuck, The Nuclear Many Body Problem 
(Springer-Verlag, New York-Berlin, 1980). 

[21 D.J. Dean and M. Hjorth-Jensen, Rev. Mod. Phys. 75 
607 (2003). 

[3] V. Zelevinsky and A. Volya, Nucl. Phys. A 731, 299 
(2004). 

[4] S. Yoshida and H. Sagawa, Phys. Rev. C 77, 054308 
(2008). 

[5] L. Ferreira, R.J. Liotta, C.H. Dasso, R.A. Broglia, and 
A. Winther, Nucl. Phys. A 426, 276 (1984). 

[6] K. Hagino and H. Sagawa, Phys. Rev. C 72, 044321 
(2005). 

[71 A. Volya, V. Zelevinsky, and B.A. Brown, Phys. Rev. C 
65, 054312 (2002). 

[8] N. Pillet, N. Sandulescu, and P. Schuck, Phys. Rev. C 
76, 024310 (2007). 

[9] N. Pillet, N. Sandulescu, P. Schuck, and J.-F. Berger, 
Phys. Rev. C 81 034307 (2010). 

[10] X. Vinas, P. Schuck, and N. Pillet, Phys. Rev. C 82, 
034314 (2010). 

[11] F.A. Janouch and R. J. Liotta, Phys. Rev. C 27, 896 
(1983). 


larger than the geometrical radius for light nuclei and 
approaches this value for heavy nuclei. Our analysis evi¬ 
denced strong shell effects. 

The pairing coherence length slowly decreases with in¬ 
creasing temperature, indicating a gradual quenching of 
pairing correlations, as is natural in finite systems. In ex¬ 
otic nuclei close to drip lines, where the Fermi energy has 
positive values, the correlation length has larger values 
and therefore the spatial correlation increases. 

The quarteting coherence length describes correlations 
between proton and neutron pairs, by overlapping their 
relative parts to the corresponding pp and nn compo¬ 
nents of the a-particle wave function. It has a similar 
behaviour, but with larger values on the nuclear sur¬ 
face. We evidenced the important role played by proton- 
neutron correlations by considering in addition the over¬ 
lap with the pn component of the a-particle wave func¬ 
tion. They change completely the behaviour of the quar¬ 
teting coherence length, namely the alpha correlation 
length has oscillating values around the a-particle geo¬ 
metrical radius. Its mean value ss 1.7 fm weakly depends 
on the nuclear mass. The analysis of the two terms en¬ 
tering the definition of the coherence length reveales the 
main difference between the pairing and quarteting cases. 
It turns out that pairing correlations are larger inside nu¬ 
cleus, while quarteting correlations are connected to the 
nuclear surface. 


ACKNOWLEDGMENTS 

This work has been supported by the project PN-II- 
ID-PCE-2011-3-0092 and NuPNET-SARFEN of the Ro¬ 
manian Ministry of Education and Research. 


[12] D.S. Delion Theory of Particle and Cluster Emission, 
(Springer-Verlag, New York-Berlin, 2010). 

[13] G. Ropke, A. Schnell, P. Schuck, and P. Nozieres, Phys. 
Rev. Lett. 80, 3177 (1998). 

[14] V. Zelevinsky, A. Volya, Phys. At. Nuclei 66 10 (2003). 

[15] A. Schiller el. al, Phys. Rev. C 63 021306R (2003). 

[16] T. Dpssing et al, Phys. Rev. Lett. 75, 1276 (1995). 

[17] G. Bertsch and H. Esbensen, Ann. Phys. (N.Y.) 209, 327 
(1991). 

[18] P.J. Borycki, J. Dobaczewski, W. Nazarewicz, and M. V. 
Stoitsov, Phys. Rev. C 73, 044319 (2006). 

[19] J. Dobaczewski, W. Nazarewicz, T.R. Werner, J.F. 
Berger, C.R. Chinn, and J. Decharge, Phys. Rev. C 53, 
2809 (1996). 

[20] M. Grasso, N. Sandulescu, N. Van Giai, and R. J. Liotta, 
Phys. Rev. C 64, 064321 (2001). 

[21] I. Hamamoto, Phys. Rev. C 73, 044317 (2006). 

[22] Id.R. Betan, G.G. Dussel, and R.J. Liotta, Phys. Rev. C 
78, 044325 (2008). 

[23] Id.R. Betan, Nucl. Phys. A 879, (2012) 14 (2012). 

[24] D.S. Delion, D. Santos, P. Schuck, Phys. Lett. B 398, 1 
(1997). 





10 


[25] D.S. Delion, M. Baldo, and U. Lombardo, Nucl. Phys. A 
593, 151 (1995). 

[26] H.J. Mang, Phys. Rev. 119, 1069 (1960). 

[27] A. Sandulescu, Nucl. Phys. A 37, 332 (1962). 

[28] D.S. Delion and R.J. Liotta, Phys. Rev. C 87, 041302(R) 
(2013). 


[29] G. Roepke, P. Schuck, Y. Funaki, et al, Phys. Rev. C 
90, 034304 (2014). 

[30] M. Mirea, arXiv:1411.3152vl (2014). 

[31] P. Moller and J.R. Nix, Nucl. Phys. A 272 502 (1995). 

[32] J. Dudek, W. Nazarewicz, and T. Werner, Nucl. Phys. 
A341, 253 (1980). 

[33] L. Liu, Z.-H. Zhang, P.-W. Zhao, arXiv:1412.5069vl 
(2014). 


